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ABSTRACT 



If binaries are common among massive stars, it will have important consequences for the 
derivation of fundamental properties like the cluster age, IMF and dynamical mass. 

Making use of the multiplexing facilities of Gemini Multi Object Spectrograph (GMOS) we 
were able to investigate the presence of binary stars within the ionising cluster of 30 Doradus. 
From a seven epochs observing campaign at Gemini South we detect a binary candidate rate 
of about 50%, which is consistent with an intrinsic 100% binary rate among massive stars. We 
find that single epoch determinations of the velocity dispersion give values around 30 kms -1 . 
After correcting the global velocity dispersion for the binary orbital motions, the "true" cluster 
velocity dispersion is 8.3 kms -1 . This value implies a virial mass of about 4.5 x 10 5 Mq or 8 
percent of the mass calculated using the single epoch value. The binary corrected virial mass 
estimate is consistent with photometric mass determinations thus suggesting that NGC 2070 is a 
firm candidate for a future globular cluster. 



Subject headings: Stars: early-type - stars: kinematics 
general - galaxies: Magellanic Clouds. 



binaries:spectroscopic - galaxies: clusters: 



1. Introduction 

30 Doradus in the LMC is the nearest available 
example of a young and massive starburst clus- 
ter. Given its proximity it is possible to perform 
a highly detailed study of its stellar component. 
The large number of massive stars present in this 
single cluster allows the statistical analysis of sev- 
eral parameters at a level of significance that is not 
available in the local and smaller Galactic clusters. 

Although a large amount of work has been de- 
voted to the study of the 30-Dor ionizing cluster, 



some important topics, like the binary fraction 
among its massive stars remain relatively unex- 
plored. If binaries are common among massive 
stars, it will have important consequences for the 
derivation of fundamental properties like the clus- 
ter age, IMF and dynamical mass. In particular, 
the orbital motions of massive binary stars in a 
young stellar cluster can have a large impact on 
the measurement of its global stellar velocity dis- 
persion and consequently on the dynamical mass 
determination. 

Bosch (1999) and collaborators pointed out 
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that the observed radial velocity dispersion of stars 
within the cluster, can be strongly affected by the 
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orbital motions of massive binaries. Single epoch 
mid r esolution spectra ob tained at ESO with the 
NTT (|Bosch et al.lll999ah were used to perform a 
radial velocity analysis on the stars that conform 
the ionising cluster of 30 Dor ( Bosch et al. I l200lh . 
The most important results of that analysis can 
be outlined as follows: 

• A high value of the velocity dispersion (a ~ 
35 kms -1 ) was obtained for the OB stars 
within NGC 2070, the ionizing cluster in 
30 Dor. 

• Estimates of the influence that an under- 
lying binary population could have on this 
determination revealed that orbital motion 
within the binary pair could mimic this large 
velocity dispersion. 

• From the photometric mass of 30-Dor a 
rather uncertain value for the virial motions 
was estimated to be around 10 kms -1 . 

Although these results are consistent with a 
100% binary population, it is necessary to confirm 
this assumption with direct evidence. The issue of 
binary frequency is still subject to debate, as most 
of the statistics o n binary sta rs are biased towards 
later type stars ( Ladal Eooih . There is evidence 
pointing towards a relatively high number of b i- 
naries among early-type stars (jMason et al . 1998). 
but previous studies based on single Galac- 
tic clusters u sing similar techniques (such as 
(jSana et al.ll2008h ) deal with smaller num- 
ber of stars, each of them observed on dif- 
ferent number of epochs. 

Here we present a new set of observations of 
NGC 2070 obtained with the Gemini Multi Object 
Spectrograph (GMOS) at Gemini South. These 
comprise multi object optical spectroscopy of 50 
early-type stars observed at least in six different 
epochs. The aim is to detect spectroscopic binary 
stars from variations in their radial velocities. 

2. Observations and Data Reduction 

Observations were performed at Gemini South 
Observatory (Proposals GS-2005B-Q-2 and GS- 
2006B-Q-21) using two multislit masks. Targets 
within each mask were selected from a previ- 
ous imaging run with GMOS according to their 



spectral types as determined and compiled in 



Bosch et all (1999a). 



Spectra were obtained during nights of September 
and December 2005 for the first run and Decem- 
ber 2006 and January 2007 for the second run. 
Approximate HJD for the listed observations are 
shown in table [TJ The instrument was set up with 
the B1200 grating R ~ 3700 centred at about 4500 
A which yields a resolution of 0.25A per pixel at 
the CCD. Although the wavelength range varies 
slightly in MOS spectroscopy, the region from 3900 
A to 5500 A is covered by our spectra. The to- 
tal integration time was split in three to allow for 
the wavelength dithering pattern needed to cover 
the gaps between GMOS CCDs. Overall signal 
to noise ratios are above 150. These ratios are 
measured on reduced spectra, dividing the aver- 
age value of the stellar continuum by the scatter 
over the same spectral range. No flux calibration 
stars were observed, as they were not needed for 
our purposes. 

Data were reduced following standard proce- 
dures, using the GMOS reduction tasks within the 
Gemini package. 

3. Radial Velocities 
3.1. Zero-point errors 

We were able to check for the presence of zero- 
point errors in our radial velocity determinations 
using the nebular lines present throughout the 
region. Nebular spectra were wavelength cali- 
brated together with the stellar spectra, as to pro- 
vide a strong template to check for variations of 
our radial velocities zero-point at different epochs. 
Nebular emission lines have strong narrow pro- 
files, enabling us to cross-correlate nebular spectra 
obtained on different nights, using fxcor within 
IRAF, that yields accurate determinations of ra- 
dial velocity differences, if present. Radial veloc- 
ities derived for each night (V ne t,) are very stable 
when observations on different epochs are com- 
pared and their variations |AV„ e f,| were found to 



X IRAF is distributed by the National Optical Astronomy 
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be negligible ( |AV4eb| < 1.0 kms \ a neb = 3.1 
kms -1 ) which suggests that there is no system- 
atic shift introducing spurious variations of stellar 
radial velocities between epochs. 



3.2. Measurements 

Radial velocities were derived measuring ab- 
sorption line profiles with the aid of the ngaussf it 
task within the STSDAS/IRAF package, follow- 
ing a simila r proc edure as the one described in 



Bosch et al 



(|200lh . This allowed us to derive in- 
dividual radial velocities for each spectral line, and 
handle each element (or each ion in the case of He) 
separately. Stellar radial velocities were derived 
using the best set of lines available, according to 
the star's spectral type. Although this means we 
are not using the same set of lines for the whole 
sample, we strictly kept the same set of lines for 
the same star on different epochs when looking 
for radial velocity variations. Uncertainties intro- 
duced when fitting individual Gaussians to the ab- 
sorption profiles are of the order of 5 kms -1 . 

Table[2]lists the complete set of radial velocities 
determined for the sample star s. Stars are lab elled 
following the nomenclature of Parkerl |l993) and 
data columns include average radial velocity (and 
its uncertainty) for each epoch. The errors listed 
in columns 3, 5, 7, 9, 11, 13 and 15 in Table [2] 
correspond to the internal error aj which is the 
quadratic sum of the standard deviation of the es- 
timated average from the set of available lines and 
the minimum uncertainty (3.1 kms -1 ) as derived 
from high S/N narrow nebular emission lines. Col- 
umn 16 shows the ratio between the standard de- 
viation of stellar radial velocity between different 
epochs and the average of 07. Spectral types are 
given for reference in column 2. 
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01 


Jan. 


2007 
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Jan. 
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Table 1: Log of observations. Observing dates are 
listed in column 1 and Heliocentric Julian dates 
for masks I and 11 are shown in columns 2 and 3 
respectively. 
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Table 2 



Stellar ide n tificat ions by IParkerI (|1993l ) are in column 1 and spectral types from 
Bosch et al.I (|1999aI ) are included in column 2. Columns 3 through 16 include radial 



VELOCITIES AND UNCERTAINTIES (ALL IN KM S 1 ) FOR EACH EPOCH OF OBSERVATION. COLUMN 17 
LISTS THE RATIO BETWEEN EPOCH-TO-EPOCH VARIATIONS AND THE AVERAGE UNCERTAINTIES. 
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4. Results 

4.1. Detection of Binary Candidates 

From the different epochs' radial velocities we 
can check for the presence of variations with time. 
To quantify this, we follow the standard proce- 
dure of comparing the dispersion of the average 
radial velocity for each star with the average un- 
certainty in the determination of each radial veloc- 
ity. This can be done calculating the "external" 
to "intern al" velocity dispe rsion ratio (cte/ci) as 



to intern al velocity dispersion ratio (<t^/(Ti) a: 
defined bv lAbt et alj (1972). Radial velocity vari 



ables can then be easily flagged out as they show 
cte/ci above 3, which is similar to say that the 
variation in radial velocity is 3 a above the ex- 
pected uncertainties. 

Binary motion is not the only source for scatter 
in the determination of radial velocities for a star. 
Stellar winds and atmospheric turbulence - both 
present in early type stars - also introduce appar- 
ent changes in the derivation of radial velocities 
at different epochs, but in these cases the inter- 
nal errors are large too. A clear example of this 
i s Star #205 in the ionizing cluster of NGC6611 
((Bosch et alJll999bl) . 

In addition to the results presented in Table [21 
the high signal to noise ratio of our spectra has also 
allowed us to detect several double-lined binaries 
that show evident variation of their absorption line 
profiles from epoch to epoch. Examples of these 
can be seen in Figures [T] and [2j Individual radial 
velocities are difficult to measure for these cases, 
therefore some of these stars are not included in 
Table These ten stars (including four present 
in Table 2) are listed in Tabled 

An inspection of Table[2]shows that 17 out of 46 
stars show radial velocity variations. If we add the 
stars present in Table [3] (excluding P613 which al- 
ready shows radial velocity variations in Table 
we have our complete set of binary star candidates, 
which rises to 25 out of 52 stars (~ 48%). 

The mass ratio of binary components, orbital 
parameters, number of observations and radial ve- 
locity accuracy, can impact the statistics of de- 
tected binaries in star clusters. In an effort to 



quantify these effects, iBosch fc Mezal (|200ll) per- 
formed Montecarlo simulations including all the 
parameters listed above. Feeding our GMOS data 
parameters into these simulations, we find that our 



detection rate is consistent with a population of 
100% binaries. 

4.2. Cluster Kinematics 

After removing the confirmed and candidate 
binaries, the non-binary population of our sam- 
ple decreases to 26 stars, still enough to calculate 
a representative value of the stellar radial veloc- 
ity dispersion (ay). In calculating it, we must 
keep in mind that each radial velocity measure- 
ment has its intrinsic error ((Tint) which must be 
subtracted quadratically from the observed (<J b s ) 
velocity dispersion. The former is the average of 
the (<7i) values over the whole sample of non ra- 
dial velocity variable stars while the latter is di- 
rectly derived as the standard deviation around 
the average radial velocity derived for the same 
subset of stars in the cluster and it therefore fol- 
lows that ay 2 = Cobs 2 — (cmt) 2 - For °ur sample, 
Cobs = 10.3 kms , (cri n t) = 6.2 kms -1 , which 
yield a value of 8.3 kms -1 for the actual radial 
velocity dispe rsion. This seems to confirm the 
suggestion by IBosch et all (|200lh based on sim- 
ulations, that the large values derived for the stel- 
lar velocity dispersion were most probably due to 
the presence of binaries. As expected, if we derive 
the radial velocity dispersion from an individual 
GMOS mask observation for a single epoch, we 
find values as high as 30 kms -1 , consistent with 
what was previously found in Bosch et al. (2001) 
from single epoch NTT data. 

This is the first time that a radial velocity dis- 
persion is found for massive stars in NGC 2070, 
free from contamination by motion of binary stars. 
We can therefore safely proceed to compare it with 
the stellar kinematics expected for 30 Dor in or- 
der to check if the cluster is virialised. Following 
th e estimations f or the photometric mass derived 
in Selman et al. (1999), we find that the stellar 



mass within the volume covered by our spectro- 
scopic survey is about 4.5 x 10 5 M Q . If the stel- 
lar velocity distribut i on is isotropic, we can follow 
Binnev &: Tremainei (1987) and estimate the ex- 
pected radial velocity dispersion as, 



o r 



4.5 x 10~ 3 M 
3Rc 



where masses are in M©, distances in pc and a 
in kms -1 . The velocity dispersion derived for 
the total mass of 30 Dor, for a half-mass radius 
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Table 2 — Continued 
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Fig. 1. — Spectra obtained for Parker 1024 on two 
different epochs, showing the variation of observed 
profiles for selected absorption lines. Both spectra 
are continuum normalised but shifted in the y-axis 
by 0.05 continuum units for comparison purposes 
and identified by their respective observing date 
during 2005. Relevant absorption profiles are la- 
belled and double troughs are identified indicating 
their blue (b) and red (r) components. Note the 
evident changes in the relative intensities of the 
blue and red troughs, particularly for Hen 4542A 
and Hen 4686A, suggesting the binary pair was ob- 
served on opposite orbital phases. 
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Hel Hel 
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Hel Hel 



4380 4400 4420 4440 4460 4480 

Fig. 2. — Same as figure [1] for Parker 260 with 
a 0.2 continuum units shift between spectra. In 
this case the profile variations show changes from 
a well-behaved absorption profile (Dec 24) to an 
asymmetrical profile (Dec 19). These changes are 
particularly evident for Hel 4387A and 4026 A pro- 
files. Although the most evident changes are seen 
as a double peaked profile in Hei447lA the pres- 
ence of nebular emission contamination in this line 
cannot be ruled out. 
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of 14 pc is about 8 krns" 1 . Our derived value of 
8.3 kms" 1 is well within the observational - and 
sampling - uncertainties. If there are more bina- 
ries yet to be discovered a slight reduction in the 
observed value is expected (although it is unlikely 
that they would show such large radial velocity 
variations that would enhance the current value 
dramatically) . 

Another interesting issue is that the single 
stars population does not confirm the evidence 
towards ma s s segr egation that was suggested in 
Bosch et al. (2001), at least for the massive stars. 
However, the fact that massive stars seem to have 
a relatively low velocity dispersion cannot rule out 
dynamical mass segregation with low mass stars. 

5. Conclusions 

Based on seven epoch observations with GMOS 
in Gemini South, we have presented observa- 
tional evidence that shows that the binary fraction 
among massive stars in NGC 2070 may be very 
high. Indeed we already detect a 50% binary can- 
didacy with only three epoch observations, which 
suggests that it is only a lower limit. The evidence 
pointing towards a high binarity fraction has an 
important effect on the massive end of the stellar 
cluster IMF, and shouldn't be overlooked. 

If a high binary fraction is common among 
massive young clusters, virial masses obtained 
from the observed global velocity dispersions 
in unresolved young clusters might be over- 
estimated by a large f actor as discussed in 
Kouwenhoven fc de Griil (|2008h . 

Regardless of the origin of the stellar radial ve- 
locity variability, the analysis of the non-variable 
subset provides an important result regarding 
kinematics of the stellar cluster itself. The radial 
velocity dispersion determined for the 30 Doradus 
ionising cluster agrees, within observational er- 
rors, with the stellar kinematics expected if the 
cluster is virialised and its total mass is derived 
from the photometric plus ionised gas masses. 
This suggests that the stellar cluster is far from 
disruption and stands as a firm candidate for a 
future globular cluster system. 

More observations in further epochs should al- 
low us to confirm the binary nature of radial ve- 
locity variables detected in this work and to derive 
orbital parameters and individual masses of mem- 



bers of the binary pairs. The success of GMOS 
in the MOS mode makes it a very efficient tool 
for discovering massive binaries in large numbers 
using a relatively small amount of telescope time. 
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hi 



Comment 



260 Asymmetrical on 19/12 

613 Hen lines look double on 24/12 

674 Double lines on 22/12 & 24/12 

805 Double lined on 23/12 

955 Looks double on 20/12 & 23/12 

956 Looks Asymm. on 23/12 
1024 Double profile on three epochs 
1031 Double lined on 07/09 & 20/12 
1035 Hei lines look broad on 24/12 
1938 Broad lines on 22/12 



Table 3: Stars with evidence of binarity f rom in- 
dividu al spectra. Stellar identifications by IParker 
() 19931 ) are shown in column 1. Column 2 describes 
the characteristics found among spectral features 
that suggest binary nature. 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 
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